Coherent ratchets in driven Bose-Einstein condensates 
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We study the response of a Bose-Einstein condensate to an unbiased periodic driving potential. 
By controlling the space and time symmetries of the driving we show how a directed current can 
be induced, producing a coherent quantum ratchet. Weak driving induces a regular behavior, and 
both space and time symmetries must both be broken to produce a current. For strong driving the 
behavior becomes chaotic and the resulting effective irreversibility means that it is unnecessary to 
explicitly break time symmetry. Spatial asymmetry alone is then sufficient to produce the ratchet 
effect, even in the absence of interactions, and although the system remains completely coherent. 
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Introduction - The physics of ratchets, systems that 
exhibit directed motion in the absence of an external 
bias, has undergone extremely rapid development in re- 
cent years [HQ. The concept is very general, and ranges 
from new technological forms of manipulating and direct- 
ing matter at nanoscale levels, to understanding how sys- 
tems in nature such as biological molecular motors func- 
tion. Fundamentally a ratchet must satisfy two essential 
requirements: the system must be driven out of equilib- 
rium by an external force, and the relevant space/time 
symmetries, which would otherwise forbid the generation 
of directed currents, must be broken. 

A well-known example is provided by a Brownian par- 
ticle in a periodic potential. Driving the system from 
equilibrium by either pulsing the potential ("flashing 
ratchet" ) or tilting it ( "rocking ratchet" ) produces a cur- 
rent if the spatial or temporal symmetry of the driving 
force is broken. In the most commonly studied over- 
damped regime, the ratchet current arises from the rec- 
tification of random fluctuations, and accordingly noise 
and dissipation are essential ingredients. This is not true 
in general, however,and surprisingly it has been shown 
recently [J 0, @, @, 0, B 0| that ratchet effects can also 
occur even in completely coherent systems. 

Considerable progress in this direction has been made 
by considering the quantum kicked rotor. Experimen- 
tally this system can be realized extremely well in gases 
of ultracold atoms held in pulsed optical lattices. While 
it was originally thought that ratchet effects would only 
arise in systems with an underlying mixed classical phase- 
space Q , recent work has shown that they can also arise 
when the phase-space is globally chaotic. In Ref. [4[ 
a quantum Hamiltonian ratchet of this type was studied 
both theoretically and experimentally, in which the cur- 
rent arose from the generation of an asymmetry in the 
momentum distribution due to the desymmetrization of 
the system's Floquet states. An alternative scheme, de- 
veloped in the context of quantum maps Q, is to use 
interference effects to produce an imbalance in the phase- 
space distribution. Quantum resonances, where the pe- 
riod of the kicks is matched to the inverse recoil velocity 



of the optical lattice, have also been proposed @, 0] as a 
means of producing ratchet accelerators. 

In this work we consider an optically trapped Bose 
Einstein condensate (BEC), since the macroscopically 
protected coherence and excellent controllability of these 
systems make them ideal subjects for investigating quan- 
tum transport effects. Instead of kicking the system 
[H E3| j we use a smoothly varying potential and so can 
expect to produce less heating effects, which we verify 
explicitly by evaluating the fraction of non-condensed 
atoms. By choosing a form for the driving which enables 
us to separately break space and time symmetries, we 
find that we can induce a directed current in a BEC start- 
ing from a symmetric initial state. This occurs by two 
distinct mechanisms; for weak driving the system under- 
goes regular oscillations, and both space and time sym- 
metries in the driving must be broken. Conversely, for 
strong driving, the system's dynamics becomes chaotic, 
and this produces an effective irreversibility which means 
it is not necessary to explicitly break the time symmetry. 

Model - We consider a BEC confined in a toroidal trap 
111 ] with a cross-section much smaller than the trap's ra- 
dius, R. The system can thus be described by an effective 
one-dimensional Gross-Pitaevskii equation (CPE) 

1 d 2 



H(t) = 



2dx 2 



+ g\^x,t)\' + KV(x,t), (1) 



where x parametrises distance around the trap, and we 
measure all energies in units of h 2 /2mR 2 . The short- 
range interactions between the atoms in the condensate 
are described by a mean-field term with strength g, and 
the condensate is driven by a time-periodic external po- 
tential with zero mean by modulating the amplitude of 
the optical potential. The archetypal form of a ratchet 
potential [2] is V = sin(x) + asin(2a; + <fi), where V is 
symmetric for <f> = tt/2, and is maximally asymmetric 
for (j> = 0, 7T. We make the unusual choice of factoring 
the potential into separate spatial and time-dependent 
components of this form V(x,t) = V(x)f(t) 



V(x) — sin(x) + a sin(2x), 
f{t) = sin(wi) +/3sin(2wi). 



(2) 
(3) 
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It is important to note that this potential does not cor- 
respond to either a purely rocking or flashing ratchet, 
and has the appealing feature of allowing the space and 
time symmetries of the system to be controlled indepen- 
dently. We plot the form of V(x) in the inset of Fig. 
Hk, and show that for non-zero values of a the potential 
becomes skewed, breaking both inversion symmetry and 
shift symmetry A static spatial potential of this type has 
been recently studied experimentally in [121 ]. 

Results - To probe the behavior of the system we eval- 
uate I(t) = (ip(x, t)\p\ip(x, t)) as a measure of the current 
flowing in the ring. As initial condition we choose the 
zero-current state ip(x, 0) = l/y/2n, which is convenient 
for experiment as it is the ground state of the undriven 
Hamiltonian, and so can be prepared by cooling. We 
numerically integrate the wavefunction in time using a 
split-operator method, and henceforth we shall take the 
driving frequency lu — 1. In Fig. [lji we show the time- 
averaged current, 7, obtained by integrating the system 
over 100 driving periods. For weak driving we can note 
a sharp peak centered at K = 0.15. As K is increased 
from this value the current drops, becomes negative, and 
passes through a negative peak at K = 2.4 This second 
peak is associated with an enhancement of (negative) cur- 
rent over a rather broad range of driving amplitudes. 



We show the time-dependence of the current for weak 
driving in Fig. [TJd. It is clear that the sharp peak in / 
is associated with regular oscillations of the current. For 
g = this oscillation is sinusoidal with a period much 
larger than T = 2tt/lo. Increasing g initially slightly 
enhances the amplitude of the oscillation, while also de- 
forming its waveform. As g is increased further, however, 
the oscillation's amplitude becomes abruptly suppressed 
(see Fig. [2ji). Examining the evolution of the system in 
detail reveals that the oscillation occurs chiefly between 
the states |0) and |2), where \n) denotes an eigenstate of 
the undriven Hamiltonian with (quantized) momentum 
nh. While a symmetric driving would equally populate 
states with positive and negative momentum, producing 
zero current, interference effects cause the asymmetric 
form to preferentially drive the system to | + 2) rather 
than | — 2), thus inducing a net current. A perturbative 
study of the Floquet states explains this result [l3| • 

To qualitatively study this phenomenon, let us trun- 
cate the wavefunction to just the two components of in- 
terest, ip{x,t) = A + Bexp(2ix), where \A\ 2 + \B\ 2 = 
l/v27T- Under the action of the Hamiltonian (Eq. [IJ, 
this yields the equation of motion for the expansion co- 
efficients x — (A B) 
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X, (4) 



where Uj are the Pauli matrices. We may now visualize 
the time-evolution of the system using the Bloch sphere 
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FIG. 1: (a) Time-averaged current (averaged over 100 peri- 
ods) a a function of driving strength K for symmetry param- 
eters a = f3 = 0.2 and nonlinearity g = 0.1 Clear peaks at 
K = 0.15 and K = 2.4 are marked with arrows. Inset: Plot of 
the driving potential V(x). For a = (black dashed curve) it 
is symmetric, but for a — 0.2 (solid red curve) it has an asym- 
metric sawtooth form, (b) Current induced for weak driving 
(K — 0.15). For g — the current exhibits a smooth sinu- 
soidal oscillation; as g is increased the oscillations are initially 
enhanced in amplitude and deviate from sinusoidal form, and 
subsequently become highly suppressed, (c) Current induced 
for strong driving (K — 2.4, g — 0.1). In contrast to the 
above case the current shows rapid quasiperiodic oscillations, 
which nonetheless maintain a stable time-averaged value (red 
dashed line). 



representation, where the north / south poles correspond 
to occupation of the states |0) / |2). For g = the Bloch 
vector will evolve under the influence of a fictitious mag- 
netic field By = aKf(t), and so will simply make a Lar- 
mor orbit in the x — z plane as shown in Fig. [2^,. This 
corresponds to the sinusoidal oscillation displayed in Fig. 
OJ,. For a larger value of g the Bloch vector will execute 
a more complicated "figure-of-eight" motion under the 
combined influence of B x and B z , shown in Fig. ^dp, pro- 
ducing the non-sinusoidal current oscillations seen in Fig. 
[T}d. When g is increased further the magnitude of the fic- 
titious field component B z is enhanced, until the Bloch 
vector is confined to making rapid circular orbits (Fig. 
Hfc) near the north pole, in a process directly analogous 
to the non-linear phenomenon termed self-trapping [ljj. 
These small orbits correspond to the low amplitude, high 
frequency oscillations seen in Fig. QJ> for large values of 
g. In Fig. [3Ji we show in detail how this transition occurs 
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FIG. 2: Top row: Bloch sphere representation of the the evo- 
lution of the effective two level model (Eq. [4} for weak-driving 
(K — 0.15). (a) g — 0: the Bloch vector makes a circular or- 
bit in the x — z plane, corresponding to a sinusoidal current 
oscillation, (b) g — 0.1: the Bloch vector makes a more com- 
plicated "figure-of-eight" orbit, producing the distortion from 
the sinusoidal waveform, (c) g — 0.15: a phenomenon analo- 
gous to self-trapping occurs, and the Bloch vector is confined 
to small orbits in the vicinity of the north pole, implying an 
interaction-induced suppression of current, (d) The suppres- 
sion of the current occurs at a very precisely defined value 
of g. The two-level approximation results (dashed red line) 
compare well with the simulation of the full GPE (solid black 
line), (e) Current produced as a function of the asymmetry 
parameters (g — 0.10). The current vanishes for a or — 0, 
and is positive (negative) in the top-left/bottom right (top- 
right /bottom left) quadrants. Contours of constant current 
are given by a/3 = cnst. in agreement with the perturbative 
result dashed lines show the current maxima. 

by plotting the amplitude of the current oscillations as 
g is increased, and demonstrate that this simple model 
indeed captures the main features of this effect. 



We now consider the case of strong driving (K = 2.4). 
In this regime the initial state is driven to a larger num- 
ber of excited levels, and consequently the two level ap- 
proximation is no longer valid. The generated current, 
shown in Fig. [TJi, instead shows an irregular, quasiperi- 
odic character corresponding to the large number of fre- 
quency components present. Despite its jagged appear- 
ance the current converges to a stable time-average after 
just a few tens of driving periods. The ratchet current 
thus arises from a completely different mechanism to be- 
fore; instead of a regular oscillation between two momen- 
tum states, the system now evolves to a superposition 
of many momentum eigenstates, which crucially has an 
asymmetric momentum distribution jij. This produces 
a dramatic difference between the symmetry dependence 
of the current in the weak and strong driving regimes. 
The symmetry properties of this form of driving were 
analyzed in Ref. [1 51 ] for a non-interacting system, and 
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FIG. 3: Ratchet current produced for strong driving (K = 
2.4) for a system with no temporal asymmetry (i.e. (3 = 0). 
Non-zero spatial asymmetry (a 0) is required to induce a 
current, and the direction of the current depends on the sign 
of a. The magnitude of the current depends only weakly on 
the nonlinearity g. 

it was concluded that both space and time symmetries 
needed to be broken to produce a ratchet current. This 
is indeed the case for weak driving, and in Fig. [5^ we 
can clearly see that the current vanishes for a, j3 = 0. In 
contrast, when the system is driven strongly only a^0 
is required, that is, it is sufficient to just break the spa- 
tial symmetry, an effect also noted in quantum resonant 
ratchets [l6[ • In Fig. [3] we show the current produced 
for strong driving when the temporal asymmetry f3 is set 
to zero. Away from a = we obtain a clear ratchet 
current whose direction depends on the sign of a. This 
occurs in analogy to the production of a ratchet current 
in a non-Hamiltonian system [TtJ , but instead of dissipa- 
tion it is the quasiperiodic evolution of the system which 
produces the effective irreversibility in time [l8j]. This 
contrasts sharply with the results obtained in Ref. [§], 
where instead the interaction was argued to play the role 
of breaking the time-reversal symmetry. As a result, the 
ratchet current we obtain has only very weak dependence 
on g (as shown in Fig. [3]), and occurs even in the non- 
interacting case. 



As well as driving the dynamics of the condensed 
atoms, the potential also has the effect of exciting atoms 
out of the BEC to form a thermal cloud. This deple- 
tion eventually leads to the destruction of the BEC, 
and seeing how rapidly this occurs allows us to assess 
the stability of the BEC under driving 10]. This can 
be done by making the Bogoliubov approximation, and 
linearizing the GPE about i ts g round-state. Following 
the Castin-Dum formalism [19], the mean number of 
non-condensed atoms at zero temperature is given by 
N(t) = J2k=io( v k(t)\vk(t)}, where (u kl v k ) are the am- 
plitudes of the Bogoliubov quasiparticle operators with 
quantized momentum k. This approximation is valid pro- 
vided that the number of particles excited from the con- 
densate is small compared to the condensate itself. 
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FIG. 4: (a) For weak driving (K — 0.15, upper black curves) 
the number of non-condensed atoms exhibits significant oscil- 
lations, arising from the oscillation of atoms in the conden- 
sate. For strong driving (K — 2.4, lower red curves) the num- 
ber rises monotonically with a similar approximately power- 
law rate of increase, (b) Occupation of non-condensed modes 
after 100 driving periods for weak driving (g = 0.1). This 
distribution mimics the spectrum of the condensate, with a 
large peak at n — +2. (c) As for (b) but for strong driv- 
ing. The distribution of momentum is again asymmetric, but 
many modes are populated. 



Conclusions - We have investigated the dynamics of 
a BEC under a periodic driving potential. In the weak 
driving regime a ratchet current can be generated by 
breaking both space and time symmetries, and inducing 
slow, regular oscillations of the condensate. For strong 
driving a more complicated quasipcriodic dynamics is in- 
duced, in which chaos acts to eliminate the long term 
memory of the system, thus mimicking thermal noise. 
Accordingly only the spatial symmetry then needs to be 
broken to produce a ratchet current. This current should 
remain stable over timescales comparable to the system's 
Ehrenfest time, which scales logarithmically with the 
number of atoms in the condensate [2(| , as also suggested 
by numerical results presented in Ref . 0] . For a conden- 
sate of 10 5 atoms we estimate this time to be of the order 
of 50 driving periods. These results expose a new vista 
of possibilities in manipulating the interplay between the 
driving potential, interactions and symmetry breaking to 
induce directed transport in quantum coherent systems. 
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In Fig. [4^ we show the time development of the non- 
condensed atoms for both weak and strong driving. In 
all cases the number of atoms does not increase expo- 
nentially, but instead follows an approximate power-law, 
N oc t 1 with 7 ~ 1.9, indicating that under the driving 
the condensate docs not exhibit dynamical instability. It 
is interesting to note that for weak driving the produc- 
tion of non-condensed atoms exhibits marked oscillations. 
These are a consequence of the large, slow current oscil- 
lations that the driving generates; the regular sloshing 
motion periodically produces peaks in the density distri- 
bution of the condensate, causing an enhanced emission 
of particles at those times and locations. As the density 
distribution can be directly imaged in experiment, this 
provides a convenient means to study the dynamics of 
the BEC in this regime. The momentum distribution of 
the non-condensed atoms also mimics the behavior of the 
condensate. In Fig. QJ) we see that weak driving essen- 
tially only excites a single Bogoliubov mode, with mo- 
mentum k = +2. As the driving is not resonant with the 
Bogoliubov frequency, however, this mode does not grow 
exponentially with time fiol ] and so dynamical instability 
is avoided. Conversely, for strong driving a larger num- 
ber of Bogoliubov modes are excited, and the momentum 
distribution is strongly asymmetric, again resembling the 
momentum spectrum of the condensate. 
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